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Abstract. A search for narrow 0(1540)^, a candidate for pentaquark baryon with positive strangeness, 
has been performed in an exclusive proton-induced reaction p + C{N) ©"""if" + C(A'') on carbon nuclei 
or quasifree nucleons at -Ebcam= 70 GeV (y^ — 11.5 GeV) studying nK'^ , pK^ and P-ft'£ decay channels of 
0(1540)''" in four different final states of the 0^ R'^ system. In order to assess the quality of the identification 
of the final states with neutron or we reconstructed yl(1520) nK% and (j) — > KiKg decays in the 
cahbration reactions p + C{N) yl(1520)7i'+ + C(A'') and p + C{N) p(j> + C{N). We found no evidence 
for narrow pentaquark peak in any of the studied final states and decay channels. Assuming that the 
production characteristics of the O'^ system are not drastically different from those of the yl(1520)-ft"'' 
and p(j} systems, we established upper limits on the cross section ratios a{0+K°) / a{A{1520)K+) < 0.02 
and a{e+K°) / a{p(f> ) < 0.15 at 90% CL and a preliminary upper limit for the forward hemisphere cross 
section a{0'^K°) < 30 nb/nucleon. 

PACS. 12.39.Mk Glueball and nonstandard multi-quark/gluon states - 13.85.Rm Limits on production 
of particles - 14.20.-c Baryons - 25.40.-h Nucleon-induced reactions 

1 Introduction seen in some experiments and not in the others was pro- 

posed in p^ . In any case it is evident that more experi- 
Although the list of the experiments supporting the first ments are needed with different beams, targets, energies 
observation m of the narrow 0+ baryon with exotic quan- and higher statistics to reject or confirm the existence of 
tum numbers is impressive [2-13], its properties (spin, par- the 6>+ baryon and estabhsh its properties. 

ity, total width and production cross sections) are not es- f^r^TTTAT^r • , , , . 

,1,.,, 1 ., -r r T,- SPHINX experiment has a long history of searches 

tablished and even its mere existence is tar from being _ . , ^ , , ? . . 

, rrii • i 1 u ^+ u lor exotic baryons and other exotic structures m van- 

proved, ihe experimental concerns about the fc*^ baryon • i i i • ■ ■ i ■ 

include the statistical significance of the observed peaks, P'^°*°';^ induced exclusive and semi-mclusive reactions. 

; ,. . J ■ J r^o As a result ot the first stage ot tfie experiment, we pub- 

evident discrepancy m its mass measured m nil ^ and pA^ .j]^ r 't f th d t' f h (M 

decay modes, the physical meaning of cuts used to enrich ^ upper imi s or e pro uc ions o ^^'^y I 

and even to see the signal. Possible sources of the false ^■^.^^^.llZ'lJ '"'rS'' ^^^^^^^^^ P^,^ 

peaks were proposed for nK+ final state [H d as weU ^^^^ ^(1385 systeins M as well as m the and ppp 

as for pK^s one M The discussion of these questions can ^^^tems [25]. Searches for the narrow Ar^(1960) baryon 

h f H |i c| |i tIi were also negative j/nlj . At the same time we found an m- 

' Ze.^^J^.eresultsine.se.rcMor OyO- started ^"of.''!; ""^^'^ '^^''^f''! .""""l 

to emerge P @, H, EH P , some of them still having ^^^^^"1^ ^he 1.7-2.1 GeV mass region. The origm of 

.1 , , T"^ r J. J.- u these structures is uncertain, they need further study and 

the status ot conference presentations. However, negative „ ... • m, r , , 

1, -iu f ii, J- J t -c confarmation m other experiments, i he status ot the stud" 

results come either from the study of specific processes „„„ u„ 1^ 1^ HTl 



/ T/>7> J v J \ e u- I, • 1 • of some of these structures can be found in 271 281 . 29l 3' 

[J/y^ and Z decays) or from high energy inclusive nu- ' — — 

cleon/nuclei interactions, whereas the positive ones come The search for the at SPHINX has its own history. 

from 'y / /v/e beams or relatively low energy proton After the publication of the paper by Diakonov, Petrov 

reactions. A possible explanation why the O'^ has been and Polyakov in 1997, the planning have started and 
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hadron calorimeter has been installed in 1998 with the 
idea to have additional capabilities to detect neutral hadror 
However, we always had in mind that in a similar model 
Weigel predicted much heavier (1580 vs 1530 MeV) 
and much wider (100 vs <15 MeV) exotic baryon. The 
relatively big width of the Weigel state could demand the 
background subtraction. After a lengthy reconstruction of 
more than 600 millions events recorded in 1998-1999, in 
2001-2002 we made a first attempt to find the 0+ baryon 
using the final states without neutral hadrons. We looked 
at the reaction^ 

p + N ^ 0+K*" + N; 0+ pK°c;, K*" ~> K-TT+ (1) 

where can be produced and compared it to the reac- 
tion 

p + N ^ pK^K*° + N; K*° ^ K+tt- (2) 

which can be used to estimate the background. Surpris- 
ingly enough, a narrow peak at M = 1548 MeV was found 
in the signal but not in the background reaction. However, 
it had a low significance (3-3.5 cr) and was found to be un- 
stable against cuts. Later on, with a better understanding 
of the neutron reconstruction, we made a quick look at 
the same reaction with different final state 

p + N + N; e+ ^ nK+, K*^* K-TT+ (3) 

and did not find expected signal. Currently we consider 
an early peak to be a normal statistical fluctuation. The 
reactions are still under study and the results will 

be available in the near future. 

In this work we present results of a search for 69+ 
baryon in more simple reaction 

p + N^e+K° + N (4) 

The basic idea of our approach is to study simultaneously 
all experimentally available final states of the Q+X'^ sys- 
tem, thus eliminating the influence of possible reflections 
and inevitable statistical fluctuations on the final judg- 
ment. 



2 Experimental apparatus 

The SPHINX spectrometer was running in the proton 
beam of the IHEP accelerator with energy Ep — 70 GeV 
and intensity / ~ (2-4) x 10^ p/spill in 1989-1999. During 
that time several modifications of the detector were made. 
The data presented in this work was obtained with the last 
completely upgraded version of the SPHINX spectrometer 
[30| . The layout of this detector is presented in Fig.^ The 
right-handed X, Y, Z coordinate system of the setup had 
Z-axis in the direction of the proton beam, vertical y-axis 
and horizontal A"-axis. The origin of the coordinate sys- 
tem was in the center of the Magnet. The main elements 
of the detector are as follows: 

^ Here and below we are using N instead of C{N) in the 
reaction notation. We do not distinguish the processes on C 
nuclei and quasifree nucleons in this work. However, they can 
be easily separated using distribution. 



1 . Detectors of the primary proton beam — scintillation 
counters Si — S4 and scintillation hodoscopes Hix^y 
and H2x,Y- 

2. The targets Ti (Cu; 2.64 g/cm^) and T2 (C; 11.3 g/cm^), 
which were exposed simultaneously. The distance be- 
tween the targets was 25 cm. The counter system around 
the target region included scintillation hodoscope 
(four counters per quadrant) and veto counters — lead- 
scintillator sandwiches Ai — A4 around the targets and 
A5 — As in the forward direction. The holes in the 
counters ^5 , Aq were matched with the acceptance of 
spectrometer. 

3. Wide-aperture magnetic spectrometer, based on the 
upgraded magnet SP-40 with uniform magnetic field in 
the volume of 100x70x150 cm^ and pr = 0.588 GeV/c, 
was equipped with proportional chambers PC, drift 
tubes DT and hodoscopes H4X , H^x , Hqx , Hj, Hsy ■ 
PC system consisted of five X and five Y planes with 
2 mm pitch and sensitive region of 76.8 x 64.0 cm^. 
DT system contained 18 planes of the thin-walled my- 
lar tubes with the diameter of 6.25 cm. There were 32 
tubes in each plane with the wires at 0°; ±7.5° with 
respect to the vertical Y-axis. The sensitivity of the 
central tube of each plane in the beam region was arti- 
ficially reduced. Electric field distribution in the drift 
tubes was described in |3^ . The space resolution of the 
DT plane was ~ 300 /xm on average. 

4. The system of Cherenkov counters for identification of 
secondary particles included RICH velocity spectrom- 
eter with photomatrix with 736 small phototubes 0i 
"3^ and hodoscopical threshold Cherenkov counter C 
with 32 optically independent channels. RICH detector 
was filled with SFg at a pressure slightly above atmo- 
spheric one. The threshold momenta of tt/ K/p in this 
detector were 3.5/12.4/23.6 GeV/c correspondingly. C 
detector with an air at atmospheric pressure had mo- 
mentum thresholds 6.0/21.3/40.1 GeV/c. Optical cells 
of C matched geometrically corresponding cells of ma- 
trix hodoscope Hi and this system can be used for 
charged particles identification at trigger level. 

5. Lead glass electromagnetic calorimeter ECAL was a 
matrix with 39 x 27 cells of 5 x 5 cm^ each. This 
calorimeter was used previously in EHS experiment at 
CERN 35]. One counter was removed for the proton 
beam to pass through the ECAL. 

6. Hadron calorimeter HCAL was a matrix of 12 x 8 
steel/scintillator 20 x 20 cm^ total absorption coun- 
ters (5Labs thickness) |3^ . 

The beam was produced diffractively off the main beam 
of the U-70 accelerator and had negligible momentum 
spread, small space dimensions (2x4mm^) and small an- 
gular divergence. This allowed an effective use of the small 
B1-B2 counters as a beam killer telescope. A special trig- 
ger logic scheme allowed the construction of up to eight 
different kinds of triggers, using as a primitive elements 
signals from scintillation and veto counters, multiplicity of 
hits in the hodoscopes and threshold Cherenkov counter, 
and total energy sum in ECAL. 
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The data acquisition system was based on MISS stan- 
dard 1^31 developed at IHEP and could record up to 4000 
events per 10 second spill of the accelerator. 

The statistics used in this work was written during 
the last run of the SPHINX experiment in March- April of 
1999. More than 600 million trigger events were written, 
corresponding to 3-10^^ live protons on targets. 



3 Data analysis 

3.1 General considerations 

The 0+ ^"-system in the reaction Q can have the final 
states presented in Tabled Typical energy of if£ for the 
processes under study is ~ 15 GeV. The decay length of 
if£ was more than 450 metres on average, thus it can 
be regarded as a stable particle and reconstructed as a 
neutral cluster in ECAL or NCAL in the same way as a 
neutron. The last two final states in Table^have only one 
charged track and are unavailable at SPHINX due to trig- 
ger restrictions, while the others were registered by three 
and five track triggers, which had the following structure: 

=ro*i?3(0-l)*H4(2-3)*il6(=3)*i/7(l-3) (5) 

T(5) -ro*if3(0-l)*i/4(4-5)*i/6(=5)*i/7(3-6) (6) 

where pretrigger Tq = S1S2S3S4 * {B1B2) * A^-s, and 
Hi{mi — 7712) means multiplicity requirement between mi 
and TO2 for number of hits in hodoscope Hi. For technical 
reasons five track trigger was implemented not from the 
very beginning of the run, resulting in a lower luminosity 
for this kind of trigger. In addition a simple beam trigger 
Tbeam, which was a greatly prescaled (« 16000) fourfold 
coincidence S1S2S3S4, was written throughout the run. 

Four experimentally available final states of the O'^K'^ 
initial state have different and to some extent complemen- 
tary properties, which are (together with the simulated 
characteristics of SPHINX detector) summarized below: 

— [nK'^]Kg: Definite strangeness, large branching ratio, 
moderate effective mass resolution, yl(1520) — > nKg 
calibration decay in the same final state, but: neutron 
should be detected, background from yl(1520). 

— [pKg]Kl: The best effective mass resolution, moderate 
branching ratio, (f> —> KgR^ calibration decay in the 
same final state, but: indefinite strangeness, Kl should 
be detected, background from (p. 

— [pKl]Kg: Moderate branching ratio, K^K^ cal- 
ibration decay in the same final state, but: indefinite 
strangeness, worst effective mass resolution, K'j^ should 
be detected, background from (f>. 

— [pK^]K^: No neutral particles in the final state (thus 
an additional constraint on the total energy in the 
event), moderate effective mass resolution, no or small 
background from the known particles, but: indefinite 
strangeness, small branching ratio and small efficiency 
(five tracks, two weak decays), lower luminosity. 



For all but [pK^]K^ final state, a neutral hadron [Kl or 
n) should be reconstructed. In the SPHINX detector ap- 
proximately 65% of neutrons and somewhat lower fraction 
of interact in ECAL, giving highly fluctuating hadron 
showers with a typical visual energy release ^ 20% of their 
total energy. All the others are totally absorbed in NCAL 
except for a small fraction (~ 5%), interacting in dead 
material in between. For neutral hadrons interacting in 
ECAL, only coordinates can be measured, the energy of 
hadron should be calculated using the exclusivity of the 
event. On the contrary, for neutral hadrons in NCAL both 
coordinates and energy can be measured in one detector, 
giving additional constraint on the exclusivity of the event. 
However, it was found by the Monte-Carlo simulation and 
verified with the study of yl(1520) nKg decays, that the 
method of neutral hadrons reconstruction in ECAL gives 
better effective mass resolution for nK~^ and pK^ systems 
due to the excellent ability of the SPHINX setup to iso- 
late exclusive processes and negligible momentum spread 
of diffractively produced proton beam. In addition, the 
sample of events with ECAL neutral hadrons is more than 
two times bigger. Only the sample of events with n{K^) 
in interacting ECAL is used in this work. The sample of 
events with n{K^) in NCAL was our strategic reserve and 
was planned to be used as a control one in the case of 0+ 
observation. 

Being products of the decay of the same initial state, 
the experimental final states (Table ^ have different final 
particle sets and even different multiplicities. In order to 
have conclusive results, a careful relative normalization 
and calibration is needed for all reactions under study, 
including the reactions selected by different triggers. It 
was provided by studying reactions 

p + N -> A{lb2{))K+ + N (7) 

and 

p + N ^p(t) + N (8) 

and then using different decay modes of 0-meson and es- 
pecially those of yl(1520)-hyperon, which has a lot of well 
measured decay modes ^8] with quite different topologies 
and multiplicities. With all this in mind, the following 
strategy was adopted: 

— STEP 1: Use beam trigger (minibias events) to in- 
vestigate the ability of the setup to isolate exclusive 
processes, the ability of the MC simulation to repro- 
duce trigger conditions for main types of trigger. High 
intensity reactions like p + N ^ piT^Tr~ + N and 
p + N ^ p7r+7r~7r° -|- N can be used for this study. 

— STEP 2: Develop Monte-Carlo generators for the sim- 
ulation of exclusive production of A{1520)K^ and pcj) 
systems (reactions (0^®) and adjust them to repro- 
duce the experimentally measured kinematics for de- 
cay modes with large statistics, yl(1520) pK^ (21k 
events) and </) K~^K~ (10k events). 

— STEP 3: Understand the ability of the setup to iso- 
late and reconstruct exclusive reactions with neutral 
hadrons {n, K^) in the final state. It can be done com- 
paring the generated and reconstructed reactions with 
yl(1520) -> nivTg and K^Kl decays. 
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HCAL 



Fig. 1. Layout of the SPHINX spectrometer: S1-S4 — beam scintillator counters (the very upstream counters Si and S2 not 
shown in this Fig.); Ti, T2 — copper and carbon targets; Ai-Ag — lead-scintillator veto counters; B1-B2 — noninterected 
primary beam veto telescope; Hix,y, H2x,y — beam hodoscopes; H3 — side hodoscope around the target; H4X, Hsx, Hgx, H7, 
Hgy — trigger hodoscopes; PC — proportional chambers; DT — drift tubes; RICH — velocity spectrometer with registration 
of the rings of Cherenkov radiation; C — multichannel threshold Cherenkov counter; ECAL — electromagnetic calorimeter; 
HCAL — hadron calorimeter (for details see text). 
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Table 1. Possible final states of the Q^if" system in the reaction It is assumed that the 0^ baryon has only two modes 
of decay with equal probabilities. Only tt+tt" decays of Kg meson are considered. 



STEP 4: Understand the ability of the setup to isolate In other words, the plan was to study the whole set of 
and reconstruct the exclusive reactions with five tracks calibration reactions 
in the final state using the decay yl(1520) yl7r+7r~. 
STEP 5: In addition, verify the ability of the setup 
to reconstruct as wide set of the different topologies 
of the final state as possible by using generated and 

reconstructed decays of /l(1520) to i^+TT^, Z'^tt" and p + N-^p[KgK^] +N (9) 
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for different purposes and in a certain order. Only after 
Steps 1-5 are successfully done, we should start to look 
for 0'^ in the set of signal reactions 



p + N 
p + N 
p + N 
p + N 



-N 
-N 
-N 
-N 



(10) 



> 

o 

d 



some of them being the part of the calibration set. If exis- 
tent and produced in the process under study, 0'^ should 
emerge simultaneously in all four final states in accordance 
with their relative probabilities and efficiencies. 



3.2 Study of calibration reactions 

The first step in the analysis is illustrated by Fig El All 
Tboam events were reconstructed with a standard track- 
ing program requiring at least one track after the magnet. 
The events with exactly two positive and one negative 
track after the magnet were then selected, and this sam- 
ple of unbiased events was used to study the efficiency of 
trigger elements and trigger itself (T(3) in this case). The 
peak in total energy at 70 GeV corresponds to exclusive 
events mainly of the pn^TT~ type and inelastic background 
after imposing trigger simulation cuts and simple addi- 
tional ones is quite small. This figure also demonstrates 
that interactions in targets can be easily isolated with a 
small background. Only events with a primary vertex in 
Carbon target were used in further study. The sample of 
events from Copper target was again left as a control one. 

To study the production characteristics of reactions iQ 
and which was the subject of the second step, the 
events of the reaction p+N [pK~^K~]+N were selected 
from T(3) trigger sample using the following criteria: 

— Two positive and one negative track after the magnet; 

— Good primary vertex in the carbon target; 

— No neutral clusters with E > 1 GeV in ECAL; 

— Energy balance, 65 < Etot = E1+E2 + E3 <75 GeV; 

— The momentum of any secondary particle > 5 GeV/c; 

— Identification in RICH as a pK^K^ system. 

The procedure of the identification of the final state by 
RICH detector was described in our previous publications 
(see Q for details) and will be not discussed here. More 
than 160000 events passed the selections cuts and results 
are presented in FigO The peaks from yl(1520) — > pK~ 
and (f) K^K~ decays can be clearly seen together with 
other well known structures. 

The events from yl(1520) and 4> peaks after background 
subtraction were used to study the production character- 
istics of the yl(1520)i^"'" and pcf) systems and to develop 
Monte-Carlo generators for the simulation of the reac- 
tions and 

The simulation of signal and normalization processes 
in the SPHINX setup was done in the framework of GEANT 
3.21 package. This included in particular: 

— Detailed description of geometry and material of the 
setup including, for example, individual wires in drift 
tubes; 
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Fig. 2. Events from the beam trigger. Upper panel: Total 
energy of secondary tracks for (l)-all two positive/one nega- 
tive track events, (2)-with trigger (7(3)) simulation cuts, (3)-no 
neutral clusters in ECAL, (4)-good primary vertex in targets. 
Lower panel: 2-coordinate of good primary vertex for case (3) 
of the upper panel. 
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Fig. 3. Effective mass distributions M{pK-) and M{K+K-) 
for the reaction p + N ^ pK^K^ + N after applying selection 
cuts described in the text. 
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— Realistic simulation of the efficiencies of trigger ele- 
ments; 

— Simulation of experimentally measured inefficiency of 
tracking devices (PC and DT) including inefficiency in 
the beam region; 

— Propagation of all secondary particles generated by 
GEANT including electromagnetic and hadron show- 
ers in the 7-spectrometer; 

— Simulation of the multibeam events with more than 
one beam particle within ±600 ns time window of the 
main process and then realistic simulation of signal 
development in various detectors. 

The simulation of the events was done up to the level of 
digitized detector responses so that they were processed 
in exactly the same way as real events. 

After a few iterations it was possible to reproduce the 
experimentally observed characteristics of yl(1520)i^+ and 
p(j) production with a reasonable accuracy. The variables 
studied were: the mass M of the system A{1520)K^ {p<P), 
the transverse momentum squared and the angles (9* , (jf) 
of the A{lb2Q) (</>) in the A{lb2Q)K+ (pcj)) rest frame (Gottfried- 
Jackson reference frame). The comparison of experimental 
and Monte-Carlo distributions (after passing the analysis 
chain) is presented in Fig. 0] for A{1520)K~^ production. 
The distributions for p(j) system are quite similar, as it 
was observed already in our earlier work j^^. Both sys- 
tems exhibit low mass enhancement and distinct evidence 
for the coherent production on the Carbon nuclei (see 
distribution) . The most significant difference between two 
systems is a cos 9* distribution, which for pcj) is more flat. 
The detailed study showed, however, that the efficiency 
for the processes under study is quite insensitive to the 
P^ (up to ~ 2 GeV^/c^) and production angles of the 
resulting NKK system. For example, for flat cos 6** dis- 
tribution for the production of yl(1520)i4r+ system, the 
overall efficiency only changes from 10.5% to 11%. 

With generators at hand, it was then possible to simu- 
late other decay channels of yl(1520) and and to compare 
the results of the simulation with experimental distribu- 
tions. 

The events corresponding to the reactions 
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Fig. 4. A comparison of reconstructed and simulated distribu- 
tions for the reaction p+Af ^ yl(1520)i(:+ -|-Af with yl(1520) 
pK~ decay. Open histogram represents data, hashed - MC. 
Only a small fraction of MC events is used in this picture. 
Also shown is the efficiency. 



N - 



■N 



(11) 



were selected from three track sample as follows: 



Good secondary decay vertex lies in the allowed region: 



-520 < 



< -270 cm. 



— Good primary vertex, composed of the vee vector, un- 
paired secondary track and beam track is in the Car- 
bon target: —555 < Zprim < —525 cm. 

— l/<7z > 3, where / = ^^scc — -Zprim and is the calculated 
accuracy for this quantity. 

— Missing energy E^miss = ^^bcam - EyT+ - E^- ~ Ek+ > 
5 GeV. 

— There is only one neutral cluster in ECAL with E > 
1 GeV. 

— The momentum of the unpaired positive track > 5 GeV/c. 

— RICH readings are consistent with the hypothesis that 
the negative and the positive tracks forming the sec- 




0.46 0.48 0.5 0.52 0.54 

M(7t+7il, GeV 



Fig. 5. Signals from decays of 7^5 tt'^tv and A — > pvr in 
the reactions p -I- A^ nKsK+ + N andp + A ^ An+TT~K+ + 
N. The fit gives: M{K^) = 496 MeV, a{K^) = 8.4 MeV, 
M{A) = 1116 MeV, a{A) = 3.8 MeV. 
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Fig. 6. Effective mass distributions M{pK~) and M{nK%) for 
tlie reaction p + N ^ [NK]K+ + TV. MC-simulated signals for 
yl(1520) are shown hatched. 
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Fig. 7. Effective mass distributions M{K^ K ) and 
M{KsKl) for the reaction p+7V ^ p[KK\ + N. MC-simuIated 
signals for (f> are shown hatched. 



ondary vertex are pions and the remaining track is 
kaon. 

— The effective mass M{tt^ti^) is within ±.2.5a of the 
Kg peak value. 

The events corresponding to the reactions 

p + N^pK^Kl + N (12) 

were selected in exactly the same way with the evident ex- 
change kaon <-> proton in the identification requirement. 
Approximately 20000 (30000) events oinKlK+ (pK^gKl) 
type were thus selected. Mass spectrum tt+tt" for the 
nKgK^ system before the final Kg selection cut is shown 
in Fig. El 

The energy of n (or K^, with the evident changes in 
formulas) is calculated as En — £^bcam — Ej^o ~ Ex+ and 
direction of flight of neutron is calculated using coordi- 
nates of the ECAL cluster and primary vertex. The ef- 
fective mass of riK'g system is then calculated using the 
tabular value of K'g mass as 

M{nKl) = M{mi+TT-) ~ Af(7r+7r") + M{KI). (13) 

In a threshold region this gives an improvement in the 
resolution by 20-30%, in agreement with the MC simula- 
tion. The resulting mass spectra M{nK^) and M {K^K'j^) 
are presented in Fig. and [7| We observe a good quali- 
tative agreement in the form of spectra for different final 
states and also between resonance signals in the data and 
Monte-Carlo. The agreement is also good numerically, as 
will be shown later. 

The reaction 

p + N ^ yl(1520)X+ + N] yl(1520) -> Att+tt- (14) 



is a main calibration process for five track trigger. It al- 
lows to cross check the simulation for three and five tracks 
events and connect pK^K^ decay mode of the O^K^ sys- 
tem to all the others. The events were selected from five 
track trigger sample as follows: 

— Five secondary tracks with YliQi — 1- 

— No neutral clusters with E > \ GeV in ECAL. 

— The total energy of all five particles corresponds to the 
energy of the incident proton: 65 < i?tot < 75 GeV. 

— Only one good secondary vertex in the allowed region. 

— Good primary vertex composed of vee tracks and three 
others tracks is in the Carbon target. 

— l/a, > 3. 

— RICH readings are consistent with the hypothesis that 
one of the unpaired positive tracks with momentum 
> 5 GeV is a kaon. 

— The effective mass of p7r^ combination is within ±2.5cr 
of the A peak value. 

Events of pK^Kg final state were selected in a similar 
way, demanding two (instead of one) different secondary 
vertices and the lone unpaired positive track identified as 
being consistent with proton hypothesis. Both tt+tt^ com- 
binations were then required to be within ±2.5(t of the Kg 
peak value. 

The quality of A signal for Att^tt^ K^ final state can 
be seen in Fig. [S] and mass spectra for some subsystems 
are shown in Fig|Sl The reaction is dominated by the pro- 
duction of the Z'(1385)*7r^i^~'" system. The peak from 
yl(1520) in the Att^tt^ effective mass spectrum is clearly 
seen and the number of events is sufficient to make quan- 
titative conclusions. Fig. [^represents the summary of the 
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Fig. 8. Effective mass distributions M{Att^), M{Atv~) and 
M(yl7r+7r") for the reaction p + N ^ ATr+n~K+ + N. The 
peaks of £'(1385)''' and 17(1385)" are clearly seen (upper pic- 
tures). The arrow shows the nominal mass of yl(1520). 



results in the study of A{1520)K~^ production in differ- 
ent decay modes. Two more decay modes are included, 
yl(1520) U+TT- and yl(1520) -> Z:°tt°, with S+ ptt° 
and U'^ —> Aj, A — > pTr~ decays, corresponding to the 
topologies "kink-l-27" and "vee-l-37" . Details of the data 
processing for these modes will not be discussed here. This 
figure can be used to estimate the systematics of the mass 
scale for quite different topologies and final particle sets. 

Up to now the results were presented at a qualita- 
tive level. Now we turn to the discussion of the numerical 
results for calibration processes. To estimate the num- 
ber of events, the mass spectra in Figs El and [3 were 
fitted by a sum of resonance and smooth background. 
The peaks of yl(1520) and were described by relativis- 
tic Breit-Wigner function with orbital momenta L = 2 
and L = 1 smeared by gaussian resolution. Widths of 
the resonances were fixed by world average values. Two 
other peaks in the pK~ mass spectrum are known to con- 
tain a lot of A*/X!* states and were described by two 
simple Breit-Wigner functions with free parameters for 
mass and width. The background was chosen to have a 
form of Pi{AMy^ ■ cxpi-P^AM - P^AM^) with free 
parameters Pi, AM = M — Afthr and threshold mass 
Afthr = Mn + Mk- The resolution parameter of gaus- 
sian function was either free or fixed by MC value. The 
results of fits are presented in Table El 

The errors in the number of events are statistical. The 
systematical uncertainties include an uncertainty associ- 
ated with the background description (polinomial func- 
tion instead of exponential with a smaller fitting range in 
mass), gaussian resolution (experimental or MC) and un- 



Fig. 9. Mass spectra for different decay modes of A(1520) in 
the reaction p + N -* A(1520)A'+ + N. The arrow shows the 
nominal mass of A(1520). 



certainties in the PDG parameters of the yl(1520) and 
resonant width. Added in quadrature, these do not exceed 
4 (8)% for the modes with high(low) statistics. 

As a first check of the consistency of the results we can 
calculate the relative branching for different decay modes. 
For yl(1520) we have 



BR[yl(1520) nK°] 



BR[yl(1520) 



2 • N[A{lb2Q) 



nKl]/BR[Kl 



pK-] 



A^[yl(1520) pK-] 
_ 2-2490/0.686 10.5 
^ 21200 "3^ 



0.94 ±0.08, (15) 



which is very close to the unity, as it should be. For 
meson decays we have 



BR[0 ^ K^Kl] 
BR[0 



N[4> ^ Kp<l]/JiJi[K''s 



N[(j} -> K+K-] 

_ 1440/0.686 12.1 
10660 ' ~iJ 



K+K-] 
= 0.64 ±0.06, (16) 



which again is close to the tabular value 0.693 ±0.018. We 
can conclude, therefore, that the processes with neutron 
and in the final state can be reproduced at SPHINX 
with the relative accuracy better than 10%. This is also 
true for other decay modes of yl(1520) not shown in the ta- 
ble, including the "five track" calibration decay 71(1520) 
A'K~^'iT~ , however with a little bit worse accuracy. It should 
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Particle 


Final state 


Events 


e, 


Resolution 


^mcas (^corr) i 








% 


cr, MeV 


nb / nucleon 


yl(1520) 


[pK-]K+ 


21200±300 


10.5 


8.1(7.1) 


1015(1400) 




[nK%]K+ 


2490±90 


3.8 


9.8(8.8) 


965 


</>(1020) 


[K+K~]p 


10660±190 


12.1 


3.4(3.5) 


202(279*) 






1440±80 


3.7 


7.8(6.6) 


188 



Table 2. Results of the fits of mass spectra in Figs. ElandQ Also the efhciencies (assuming YiYi{K% ■k'^-k~) — 100% for the 
decay modes with K%) and experimental(Monte-Carlo) resolution are shown. The number of events corresponds to fits with 
experimental resolution. For the meaning of cross section values see text below. 



be noted, in addition, that Monte-Carlo and experimental 
effective mass resolutions are close to each other, the MC 
resolution being typically a little bit better. This is impor- 
tant in a search for narrow states, where the estimations 
rely heavily on the resolution. 



Another question of major importance is the abso- 
lute calibration of cross sections. In most of the models 
the cross section for the production of the baryon 
in a wide class of exclusive processes is directly propor- 
tional to it's width, being governed by the same constant 
90NK- Cross section for A{1520)K~^ and pcj) production 
was calculated as cr = {1/L) ■ N/ (BR- e), where luminosity 
L per nucleon was estimated assuming A^^^ dependence 
of the cross section on the mass number. The luminos- 
ity for Carbon target was found to be 884 events/nb for 
r(3) and 445 events/nb for T(5y The measured cross sec- 
tions (Tinoas Can bc found in Tabled The data is prelimi- 
nary, as not all corrections, common to all processes and 
not influencing the relative quantities like ((T^ and 
were included into efficiency calculations. These correc- 
tions include, for example, rate dependent accidentals in 
veto counters. In general they are small and are under 
study now. Can we check the cross section values from in- 
dependent measurements? As far as we know, there are 
no data on exclusive production of A{1520)K^ system 
in nucleon-nucleon(nucleus) interactions and our result 
seems to be the first one of this kind. The existing data 
for ^-meson production in the reaction pp (ppp is scarce 
and do not allow the extrapolation to our energy. How- 
ever, we can use the data for reaction pp — > Lupp, which 
exist at higher as well as lower energies. The analysis, 
done in jSfli] . allows to extrapolate to our energy, giving 
the value (36 ± 4)/2 = 18 ± 2 /xb/nucleon for cup forward 
hemisphere production cross section. Using then cross sec- 
tion ratio a^/a^i = (1.55 ± 0.31) x 10~^, measured in our 
study of OZI-rule jifll, we arrived at the prediction for 
(pp production cross section cr^p — 279 ± 64 nb/nucleon. 
This value is close, but somewhat higher than ameas — 
202 nb/nucleon (Table|21l. Being conservative, we used this 
(higher) cross section as input. It is denoted by (*) sign 
in Tabled The correction factor fccon- = 279/202 = 1.38 
was used then to calculate the corrected cross section for 
yl(1520)i4r+ production, in Table El The same conserva- 
tive correction factor will be used later in the estimations 
of absolute cross section for production. 
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Fig. 10. Effective mass spectra of p/('£ and pK% systems in the 
reaction p + N ^ pK%Kl with {M{K%Kl) > 1.04 GeV, 
lower row) and without "0— cut" (upper row). MC-simulated 
signals (hatched) correspond to cTg+j^o /o"yi(i52o)K+ — O-l- 



3.3 Study of signal reactions 

The calibration reactions H12I) and (|ll|l are at the same 
time the reactions, where the baryon can be searched 
for. Effective mass spectra M{pKl) and M{pKg) for the 
reaction H12I) are shown in Fig. ^| and mass spectrum 
M{nK^) for the reaction (|ll|l in Fig.^] No evident struc- 
tures can be seen in these distributions, except for a hint 
on a shoulder in the pK^ mass spectrum at a mass of 
«1510 MeV. However, this structure is completely absent 
in M{pKl) mass spectrum as well as in M{nK^) mass 
spectrum. Thus we don't see the signals for the 0+ baryon 
and only upper limits can be produced from these distri- 
butions. 

We are searching for narrow signals and will assume 
that effective mass distributions are completely dominated 
by experimental resolutions for each of the systems under 
study. This assumption is valid for the width (FWHM) 
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M(nK*), GeV M(nK*), GeV 

Fig. 11. Effective mass spectrum of nK^ in tiie reaction p + 
iV ^ nK+K% + witii {M{nK%) > 1.55 MeV, lower row) 
and witiiout "yl(1520)— cut" (upper row). Pictures on tlie right 
panel are the sums of real data and MC-simulated signals for 
o"0+ j^o/o"yi(i520)if+ = O-l- Curves in the right column are the 
results of the fit of left column distributions. 



M(pKs|), GeV M(pKs,), GeV 

Fig. 12. Distributions for the reaction p + N ^ pK%K% + 
N. Upper row — M {1x^1x2) vs M{-k^'k^) with all cuts but 
K% selection cut. Lower row — 

M{pK%2) vs M{pK%{) after 
all cuts. Left column — data, right column — Monte-Carlo 
simulation of reaction p+ A ^ 0+K°+N, 0+ K° [pK'^]K'^. 



of the 0+ baryon < 6 - 10 MeV. The width like that 
and lower is indicated in most precise "positive evidence" 
experiments and the assumption doesn't seem to be hm- 
iting. To achieve the hmits on the number of events and 
then on production cross sections, the resolution and effi- 
ciency for each of the final state is needed. These can be 
estimated only by Monte-Carlo simulation. However, the 
production mechanism for 0^ baryon is unknown, so some 
model is needed to simulate its production. We simulate 
the production of O'^K'^ system as being similar to that 
of A{1520)K+ and pc/) systems (P| and M{e+K°)) and 
the distributions of the decay angles of the 0^K^ system 
in Gottfried- Jackson frame were assumed to be isotropic. 
The assumptions about production characteristics are not 
as limiting as they seem. The efficiencies demonstrate very 
weak dependence on and decay angles (see the discus- 
sion about the simulation of the A{1520)K^ production 
above) and only if prefers to originate from the O'^K^ 
system with very big mass (> 3.5 GeV) and very high Pt, 
the efficiency for the O'^K^ system detection would be 
significantly overestimated. 

The results of the simulation of an interesting case of 
the pK^K^ final state in comparison with the real data 
are shown in Fig. 1121 and one-dimensional plots, used to 
set upper limits, in Fig. 

We also made a search for the ©"''"''(1540) baryon in 
the reaction p^N ^ 0++K- + N; 0++ pK+ . This 
state emerges as a partner of 6*+ in some theoretical mod- 
els. The study was a byproduct of the investigation of the 



calibration reaction p + N ^ pK^K^ + N . The pK^ 
effective mass distributions are shown in Fig. ^1 for all 
events, for the events with "yl(1520)-cut" (the events with 
1.50 < M{pK~) < 1.55 GeV were excluded), and events 
with combined "yl(1520) -I- 0-cut", where in addition the 
events with 1.01 < M{K+K-) < 1.03 GeV were ex- 
cluded. Neither distribution shows statistically significant 
signal and only the upper limits can be set for O^^K^ 
production. 

To get the limit on the number of events from the 
distributions Figs. ll0lllllll3l and ll4l we tried two different 
methods. In the first method the distributions were fitted 
by a smooth curve plus a gaussian function with fixed 
mass and a resolution, determined by MC simulation for 
a particular final state. The 90% confidence level upper 
limits for the number of events were then estimated in a 
usual way. In the second method the events in the mass 
window near the presumed signal region were excluded 
from the fit and 90% confidence level upper limits were 
calculated as iV = max(0, Z\n) -I- 1.28Vb with An being 
the excess over estimate background b. The mass window 
was from 20 MeV for [pK'i,]Kl system to 40 MeV for 
[pKl]Kg one. To understand systematics, we varied the 
fitting range and background function. In the first method 
we also made fits with MC resolution enlarged by 1 MeV. 
Both methods give similar results with no more than 20% 
difference in upper limits. In the rest of the paper we are 
using the results from the first method. 
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Fig. 13. Effective mass spectrum of pKg system for the re- 
action p + N pKgKg + N. Lower picture corresponds to 
bin shift by 5 MeV. The MC signal corresponds to the cross 
section ratio crQ+^o/cryi(i52o)K+ = 0.1 




1.65 1.7 1.75 1.8 

M(pK"'), GeV 

Fig. 14. Search for 0'^^ baryon in the reaction p + N 
0++K-+N; 6»++ pK+. The histogramms are with no cuts, 
with yl(1520) and yl(1520)-|-(/f) cuts. The curves are results of the 
fit described in the test. The MC signal from narrow 6'^''' cor- 
responds to the cross section ratio o"g)++^- /(Tyi(]^52o)if+ ~ 0.01 
and BR{0++ pK+) = 100%. 



Our results in a search for narrow exotic baryons are pre- 
sented in Table 131 where efficiencies, gaussian resolutions, 
fitted number of events and upper limits on cross section 
for the reactions p + N ^ O^K^ + N and p + N ^ 
0++K^ + N can be found. 

The upper limits for production cross sections are pre- 
liminary. They were calculated from the upper limits on 
the number of events using the procedure described above 
for the calibration reactions. The measured values were 
corrected with the same correction factor fccorr = 1-38. 
The systematic errors for the production cross sections 
are estimated to be < 25%. 

Since different "positive" experiments reported differ- 
ent masses for the 0^ baryon, we made a scan in the 
effective mass for each of the different final states. The 
results of the scan are presented in Fig. ^1 A poor up- 
per limit for the [pi^gjifg final state in the mass range 
1515-1530 MeV is a consequence of a bump in pK^ mass 
spectrum, which is partly supported by a shoulder in pKg 
mass spectrum of [pKg]Kl final state. However, an inter- 
pretation of this bump as a (relatively wide) 0'^ baryon is 
absolutely excluded by nK'^ data. It is also possible that 
one-star Z'(1480)-hyperon really exists, producing the ir- 
regularities in pKg mass spectrum (but why not in pK^l). 

We tried to apply additional cuts in order to find a 
signal. These cuts included, in particular: 

— Cut on transverse momentum of the NKK system (co- 
herent region P|. < 0.1 GeV^/c^, high region, etc). 

— Cut on transverse momentum of the NK system. 

— Cut on the effective mass of the NKK system. 

— Cuts on the decay angles (cos6'*,0*) of the [NK]K 
system. 

None of these cuts allowed us to see a statistically signifi- 
cant signal in the region of the 0"*" baryon. In addition we 
don't see any evidence for the existence of narrow struc- 
tures in NKK effective mass distributions, in particular 
the state with M « 2.4 GeV, indicated by CLAS data 0. 
It should be noted, that with the extreme values of cuts, 
when statistics is started to be low, irregularities in differ- 
ent mass spectra begin to appear, sometimes at a mass of 
the 0^ baryon. In our case, however, they never occured 
simultaneously in the same place under the same cuts in 
more then one mass spectrum. 

Our distributions for "allowed" {S = —1), "exotic" 
(5 — +1) and "mixt" (unknown) strangeness are com- 
piled in Fig. ^1 The compilation includes pK^, nKg, 
pK^, nK~^, pKl and pKg mass spectra from pK^K~, 
nKgK~^ and pKgK^ final states. In the distributions for 
"exotic" and "mixt" strangeness events, corresponding to 
the production of 0, yl(1520), or both, were excluded. 

The relative yield 0^/71(1520) is the most commonly 
used variable to compare different experiments. In our case 
it transforms to the cross section ratio a{0^K'^)/a{A{152O)K^). 
The best way to calculate this ratio is to use the nK^K^ 
final state, which is common to both systems. Many fac- 
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Particle 


Final state 


Events 


% 


Resolution 
(a), MeV 


aeorr, 90% CL 

nb/nucleon 


(Q+n 54-fl"! 




55 ± 43 
10 ± 33 
(yli52o cut) 


3.3 
2.2 


9.8 
10.1 


< 32 

< 26 






48 ±29 
26 ±25 
((j) cut) 


3.0 
2.7 


5.8 
5.5 


< 53 

< 42 






6 ±43 
-14 ±37 
{4> cut) 


2.6 
2.3 


11.8 
11.3 


< 54 

< 39 






-4± 7 


0.9 


7.5 


< 52 


e++(1540) 


[pK+]K- 


-57 ± 100 


10.0 


8.0 


< 2 



Table 3. Upper limits for reactions p + N ^ 0+ K<^ ± TV and p + N 9++ K~ + iV at Bp = 70 GeV. The number of 
events corresponds to the fit of the distributions with smooth background function plus gaussian with M = 1540 MeV and MC 
resolution. The efficiencies was calculated assuming BR(/irs tt+tt") = 100%. 
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Fig. 15. Cross section limits on the baryon production for 
different masses of 0^. 



tors cancel in the ratio and we arrive to {A* = yl(1520)) 

cr(6)+/?0) 



Ra* 



Ne BR{A* ^ NK) e^. 
iV^ ■ BR(6)+ ^ NK) ' 



0.45 



cr{A*K+) 



Na* 80+ 



— - (17) 



where we used PDG value BR(^* NK) = 45% and 
BR(6»+ NK) was assumed to be 100%. Number of 
events and efficiencies for corresponding decay modes can 
be read from Tables |3 and [SJ and we finally get 

i?... 0.45. ^^.^=0.011 ±0.007, (18) 
2490 ± 90 0.033 ^ ' 
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Fig. 16. A summary plot with M(A''if/J^)-distributions for 
S = — 1 (upper row), 5* = ±1 (middle row) and mixt S (lower 
row). The distributions for S = ±1 are superimposed on the 
distributions for S = — 1 as hatched histograms. 



with an upper hmit Ra- < 0.02 at 90% CL for M{0+) = 
1540 MeV. The mass dependence of this limit can be un- 
derstood from [nK~^]Kg curve in Fig. El This estimate is 
free of any correction factors, and it is hard to imagine that 
the efficiency of 0^ detection is significantly lower than 
that of yl(1520). Note that for a{0+K°) « cr(yl(1520)ii:+) 
and £0+ « ea* we should see as much as 5000 6>+ nK^ 
decays. 

The ratio R^ = a{0^ K'^) / a{p(t)) can be estimated in 
the same way as Ra* and the result is R^ < 0.15 at 90% 
CL. 
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5 Discussion 

Our result for Ra* is based on the sample of ^2500 yl(1520) 
nK^ decays or, equivalently, on sa21000 yl(1520) pK^ 
decays, corresponding to « 900k of yl(1520) produced 
in the reaction p + N ^ 71(1520)/^+ + N. With some- 
what smaller number of yl(1520) HERA-B 19] reports 
(preliminary) Ra' < 0.02 for inclusive pA interactions at 
920 GeV/c and mid-rapidity («6000 yl(1520) ^ pR- de- 
cays) and ALEPH |21,| gives Ra* < 0.1 for inclusive decays 
of Z, with a sample of 2819 yl(1520) pK' decays. 

Low values of Ra* in these "negative" experiments 
are in a striking contrast with this ratio (when avail- 
able) for "positive" ones. In quasi-real photoproduction 
7*d ^ e+X HERMES found Ra- = 1.6-3.5 with « 
800 yl(1520) pR- events and SAPHIR H in exclusive 
reaction jp nR^Rg estimated a{0^Rg) w 300 nb, 
corresponding to Ra* ~ 1/3. 

A small value of Ra* , found by HERA-B and SPHINX 
in proton-induced reactions, corresponds, however, to quite 
different physical processes and kinematical regimes. It 
is natural to assume that this value should hold for any 
proton(nucleon)-induced reactions. It would be interest- 
ing to have this ratio for the "positive" experiments for 
nucleon-induced reactions. Unfortunately, yl(1520) produc- 
tion was not studied by SVD-2 collaboration |^, where 
some signal for 0"'" pKg was found in inclusive pA 
interactions at 70 GeV/c. The full cross section of 
production in the fragmentation region {XF(pRg) > 0) 
was found to be (30-120) /ib/nucleon. We can compare 
this rather big cross section with the measurements of 
the inclusive production of yl(1520) in neutron/nucleus 
interactions at 40 GeV/c where the cross section of 
yl(1520) production was found to be (70-90) ^b/nucleon 
for Xp{A{152Q)) > 0. Assuming the ratio for the inclu- 
sive cross sections of 71(1520) production for the beam 
momenta 70 and 40 GeV/c to be 1.5 (for well studied A 
production the corresponding factor is ^1.4), we conser- 
vatively find Ra* > 30/(90 • 1.5) = 0.22. 

In the other experiment 11], a positive signal was 
found in the pRg system in p -f CsHg collisions at 10 
GeV/c. Preliminary cross section for CaHg molecule was 
found to be 90 ^b, corresponding to ~7-9 ^b/nucleon (our 
estimate). The production of yl(1520) was not reported. 
We can estimate Ra* for this experiment using the re- 
sults of 42], where the inclusive cross section for yl(1520) 
production was measured innp interactions at 4-8 GeV/c 
and found to be ~ 7 /xb/nucleon. A conservative estimate 
gives Ra- > 0.2-0.3. 

Thus the results 0, ^] seem to be in direct contradic- 
tion with HERA-B and SPHINX data. 

Our upper limit on the absolute cross section produc- 
tion of 6*^ can be used to set upper limit on its width, 
though in a model dependent way. The forward hemi- 
sphere cross section for the reaction pp — > O^R^p was 
calculated in 43] using hadronic lagrangian with empiri- 
cal coupling constants and form factors. The cross section 
was found to have a maximum of 13 at ~ 4.5 GeV. 
The last point in their figure at ^/s — 7 GeV gives the 
value w6.5 ^b. We can extrapolate this value for a{0'^R^) 



to our energy {^/s = 11.5 GcV) using (conservatively) 
V-Pbcam dependence and find a{0+ R°) « 900 nb. As 
—both cross section and total width of 0'^ are propor- 
tional to g0]\jx and the width of 0~^ was assumed to be 
/^[Q+j =20 MeV in the calculations, our upper hmit on 
the production cross section transforms to the upper limit 
on the width, r[0+] < 20 • (32/900) = 0.7 MeV (the limit 
from nR^ decay mode was used). This result is in a good 
agreement with constraints on /^[^^J deduced from R^N 
and R'^d data ^^l- It should be noted, in addition, that 
with r[0+] < 0.7 MeV the calculations 43] would give 
the cross section for the reaction pp 0^S^ < few tens 
of nb, in a disagreement with COSY-TOF measure- 
ments of 400 ± 100 ± 100 nb/nucleon for this reaction at 
P=2.95 GeV/c. 



6 Conclusion 

In a high statistics experiment with the SPHINX facility 
at IHEP accelerator we have searched for the production 
of the 0+ baryon in exclusive reaction p-fA^ — > 0^R'^+N 
at the energy of 70 GeV. For the first time the search was 
done simultaneously in all possible decay modes of 0+ . We 
did not see statistically significant signals and found that 
the production of the 0~^R'^ system is very small (if any) 
compared to the production of A{1520)R~^ system and 
small (if any) compared to OZI-suppressed production of 
the p(f> system. 

The work is underway in the search for 0^ in the 
reaction p + N ^ 0^R^tt^ + N"^ . In a search for the 
exotic mechanisms for 0^ production we also plan to 
study the final states with baryon-antibaryon pair, like 
p + N ^ 0+X;+p + N. 

It is a pleasure to express our gratitude to the staff of IHEP 
accelerator and especially A. A. Aseev for the invaluable help 
with the beam adjustment in the critical moments of data tak- 
ing, and to A.M. Zaitsev for useful discussions. This work was 
partly supported by Russian Foundation for Basic Researches 
(grants 99-02-18252 and 02-02-16086). 
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